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Abstract. We report the results of two XMM-Newton observations of the ultra-compact low-mass X-ray binary 
4U 1850—087 located in the galactic globular cluster NGC 6712. A broad emission feature at 0.7 keV was detected 
in an earlier ASCA observation and explained as the result of an unusual Ne/0 abundance ratio in the absorbing 
material local to the source. We find no evidence for this feature and derive Ne/O ratios in the range 0.14-0.21, 
consistent with that of the interstellar medium. During the second observation, when the source was ~10% more 
luminous, there is some evidence for a slightly higher Ne/O ratio and additional absorption. Changes in the Ne/O 
abundance ratio have been detected from another ultra-compact binary, 4U 1543—624. We propose that these 
changes result from an X-ray induced wind which is evaporated from an O and Ne rich degenerate donor. As the 
source X-ray intensity increases so does the amount of evaporation and hence the column densities and abundance 
ratio of Ne and O. 
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1. Introduction 

ASCA and BeppoSAX observations of ultra-compact 
(Porf)<l hr) low-mass X-ray binaries (LMXBs) have 
revealed two possible spectral differences compared to 
the longer period systems. These are (1) the presence 
of a discrete spectral feature near 0.7 keV (attributed 
to Ne) in the ASCA spectra (Juett et al. [Wn\i of 
X 1850-087, X 1543-624, X0614-f01 and X 0918-549 
(the orbital periods of the last 2 sources are unknown, 
but their optical faintness is consistent with an ultra- 
compact nature) and (2) their best-fit parameter values 
when fitted with a disk-blackbody and Comptonized con- 
tinuum (Sidoli et al. 12001(1 . In the case of the ultra- 
compact sources X0512-401, X 1820-303, X 1850-087 
and X 1832—330, fits to the BeppoSAX spectra give signif- 
icantly lower disk-blackbody temperatures than for other 
LMXB. In addition, the Comptonization seed photon tem- 
peratures appear consistent with those of the inner disk 
regions. Also the LMXBs located in the globular clusters 
NGC 6652 (Parmar et al. l^niIT|l and Terzan 5 (Heinke et 
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al. IMini show very similar spectral properties, suggest- 
ing an ultra-compact nature for these binary systems (see 
Verbunt i2005< for a review) . 

The X-ray burst source X 1850-087 (Swank et al. 
1976) is an X-ray binary located in the galactic glob- 
ular cluster NGC 6712, the least concentrated amongst 
those that host a luminous X-ray source. A short period 
(20.6 minutes) UV modulation was discovered with HST 
from the Hkely optical counterpart (Anderson et al. . 1993,1 . 
implying a degenerate companion of O.O4M0 (Homer et 
al. ll996|l . The source is located '^6" or O.liO.l core radii 
from the cluster center (Hertz & Grindlay 1983). Another 
UV-excess star was discovered in the core of NGC 6712 
with the ESQ Very Large Telescope (Ferraro et al. 2000), 
a few arcsec away from the LMXB. The presence of these 
two interacting binaries inside the core of the low-density 
cluster NGC 6712 suggests that the interaction of the clus- 
ter with the disk and bulge of our Galaxy during numerous 
orbital passages plays a role in the formation of LMXBs 
in globular clusters (Ferraro et al. 2000). Moreover, there 
is evidence supporting a scenario where NGC 6712 was 
much more massive in the past and that it experienced a 
significant mass evaporation produced by the tidal force 
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due to interactions with our Galaxy (see, e.g., Paltrinieri 
et al. 120011 and references therein) . 

EXOSAT observations of 4U 1850-087 revealed a com- 
plex spectrum. The best-fit was obtained with a model 
consisting of a power-law with a photon index, a, of 
0.4 with an exponential cut-off at ~1 keV, together 
with a blackbody with a temperature, kT, of 2.4 keV 
and absorption, A'h, of <5xl0^^ cm~^ (Parmar et al. 
1989). A thermal bremsstrahlung {kT — 1.7 keV) ab- 
sorbed by 5 X 10^^ cm~^ is a good approximation to 
the ROSAT Position Sensitive Proportional Counter spec- 
trum (Verbunt et al. 1995). During a BeppoSAX survey 
of the bright LMRXB located in galactic globular clusters 
(Sidoli et al. l20t)l|l the 0.3-50 keV spectrum was fit with a 
disk-blackbody and Comptonized continuum with iVu = 
3.9 X 10^^ cm~^, an inner disk temperature, fcTin, of 0.6 
keV, an inner projected radius of ~5 km (for an assumed 
NGC6712 distance of 6.8 kpc, Harris [T996(l . a tempera- 
ture, kTo, of the input "seed" photons, of 0.8 keV (consis- 
tent with the inner disk-blackbody temperature), an elec- 
tron temperature, kTe, of 70 keV, and an optical depth, r, 
of 1.7. The 0.1-100 keV luminosity was 1.9x10^6 erg s"^ 

Analysis of the ASCA Solid-state Imaging 
Spectrometer (SIS) data (Juett et al. I2001|l from 
4U 1850—087 revealed the presence of a spectral feature 
near 0.7 keV. A good fit to these data was found with 
an absorbed kT = 0.4 keV blackbody together with a 
power-law with a = 2.1 when the relative abundances 
of O and Ne were allowed to vary. Both components are 
absorbed by iVn — 2.9x10^^ cm~^, with a relative (to 
solar) 0/H abundance of 0.37±0.06 and Ne/H abundance 
of 1.9±0.3. The authors interpreted this excess absorption 
as due to neutral Ne-rich material local to the binary. 
Preliminary results from the 0.4-2 keV 4U 1850-087 
XMM-Newton RGS spectra were reported in Sidoli et 
al. (200'P, who found no evidence for an anomalous Ne/0 
abundance ratio. Recently, analysis of Chandra Low- 
Energy Transmission Grating Spectrometer (LETGS) 
data confirmed this result (Juett & Chakrabarty |2005) , 
measuring a Ne/0 ratio of 0.22±0.05 consistent with 
that expected from the interstellar medium (ISM) of 0.18 
(Wilms et al. 2000). 

Here we report the results of XMM-Newton observa- 
tions performed in order to investigate the nature of the 
0.7 keV feature found with ASCA. We use the following 
updated values for the globular cluster NGC6712 param- 
eters (Paltrinieri et al. l2001|l : a distance of 8±1 kpc (note 
that a previous estimate for the distance was 6.8 kpc, 
Harris UM^ and a reddening E{B - V) ^ 0.33 ± 0.05. 
Adopting the relation Ay = 3.1 E{B — V), and Ay = 
Nu X 5.59 X 10^22 ^^-2 (predehl & Schmitt ITM5|l the 
optical reddening translates into an ISM column density 
of (1.8±0.2)xl02i cm-2 to NGC6712. 

2. Observations 

The XMM-Newton Observatory (Jansen et al. IMIHl in- 
cludes three 1500 cm^ X-ray telescopes each with an 



European Photon Imaging Camera (EPIC) at the focus. 
Two of the EPIC imaging spectrometers use MOS CCDs 
(Turner et al. I2001|l and one uses a pn CCD (Striider 
et al. I2001|l . Behind two of the telescopes there are 
Refiection Grating Spectrometers (RGS, 0.35-2 keV; den 
Herder et aL^J). XMM-Newton observed 4U 1850-087 
twice, due to visibility problems, in 2003 September and 
October, about 12 days apart (see Tabled for the obser- 
vation details). 

Data were reprocessed using version 6.1 of the Science 
Analysis Software (SAS). Known hot, or flickering, pixels 
and electronic noise were rejected using the SAS. The lat- 
est response matrices were used (updated to 2004-12-03, 
which should improve the agreement between the MOS 
and pn below 1 keV, Saxton i2004l . while the ancillary 
response files were generated using the SAS task arfgen. 
Spectra were selected from single events only (pattern 0) 
for the MOSl timing mode (only pattern has been cali- 
brated in this instrument mode) while for M0S2 patterns 
from to 12 and for the pn patterns from to 4 were se- 
lected. Source counts were extracted from circular regions 
of 40" radius centered on 4U 1850—087 for the pn and the 
M0S2. With the SAS task epatplot we verified that pn 
Small Window data are not significantly affected by pile- 
up, whereas in both MOS spectra pile-up was evident. 
Thus, we minimized the effects of pile-up by extracting 
M0S2 events in an annulus outside of a 10" radius core of 
the 4U 1850—087 point spread function, and MOSl events 
from a wide column outside the central 15". A comparison 
between MOSl and M0S2 spectra revealed that after this 
selection the source spectral shapes observed by the two 
instruments were similar. We use the pn for the determina- 
tion of the source fiux. Background counts were obtained 
from similar regions offset from the source position. The 
backgrounds do not show any evidence for flaring activ- 
ity, so the entire nominal exposure times were considered. 
For both observations, the RGS spectra were analyzed as 
produced by the pipeline processing performed by XMM- 
Newton Survey Science Centre. 

In order to ensure applicability of the statistic, 
the extracted spectra were rebinned such that at least 20 
counts per bin were present and such that the energy reso- 
lution was not over-sampled by more than a factor 3. Note 
that no systematic uncertainties were added to the spec- 
tra. All spectral uncertainties and upper-limits are given 
at 90% confidence for one interesting parameter. 

3. Results 



3.1. Lightcurves 

Lightcurves for the two observations in soft (0.3-2 keV) 
and hard (2-10 keV) energy ranges were extracted in or- 
der to search for variability and hardness ratio variations 
(Fig. [SJ. The source was slightly harder and more intense 
during the second observation. Thus, the two observations 
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Table 1. XMM-Newton on-axis observation log of 4U 1850-087. Two observations were performed. The MOSl, M0S2 
and pn cameras all used the medium thickness filter. 



Obs. 


Start time 




End time 


Inst. 


Net Exp. 


Mode 




(dy men yr hr:mn) 


(dy men yr hr:mn) 




(ks) 




1 


27 Sep 2003 09:05 


27 Sep 2003 12:22 


MOSl 


11.7 


Timing Fast Uncompressed 


1 


27 Sep 2003 09:05 


27 Sep 2003 12:22 


MOS2 


11.5 


Imaging Prime Partial Window 


1 


27 Sep 2003 09:10 


27 Sep 2003 12:23 


pn 


8.1 


Imaging Small Window 


1 


27 Sep 2003 09:04 


27 Sep 2003 12:28 


RGSl 


12.1 


Spectroscopy 


1 
1 


27 Sep 2003 09:04 


Z/ oGp ZUUo iZ.Zo 




19 1 


Spectroscopy 


2 


09 Oct 2003 08 


15 


09 Oct 2003 10:35 


MOSl 


8.3 


Timing Fast Uncompressed 


2 


09 Oct 2003 08 


15 


09 Oct 2003 10:39 


MOS2 


8.4 


Imaging Prime Partial Window 


2 


09 Oct 2003 08 


20 


09 Oct 2003 10:41 


pn 


5.9 


Imaging Small Window 


2 


09 Oct 2003 08 


14 


09 Oct 2003 10:43 


RGSl 


8.8 


Spectroscopy 


2 


09 Oct 2003 08 


14 


09 Oct 2003 10:43 


RGS2 


8.8 


Spectroscopy 




10 11 12 8.5 9 9.5 10 10.5 

Time (h) Time (h) 

Fig. 1. EPIC pn lightcurves in two energy bands (S, 0.3-2 keV and H, 2-10 keV) and hardness ratios, H/S, during 
the first (left) and second observation (right). Time is in hours of 2003 September 27 (for the first observation) and of 
the 2003 October 9 (for the second observation). The binning is 64 s. 



were analyzed separately. Since within each individual ob- 
servation the source does not show evidence for intensity 
or hardness variations (see Fig.^, we can safely consider 
two separate spectra extracted from each of the two ob- 
servations without making any further selections. The sta- 
tistical quality of the data, combined with the length of 
the observations, does not allow for a meaningful search 
for periods around the optical period. 

3.2. Spectra 

We performed separate spectral analysis for the two 
XMM-Newton observations. We first studied the pn, 
MOSl and M0S2 spectra in the energy range 0.3-12 keV. 
We noticed that the pn spectrum, in both the observa- 
tions, showed a significant departure from the MOSl and 
M0S2 shapes, especially around 0.6 keV, where a large 
excess is present only in the residuals of the pn spectrum 
(Fig- 13 • A significant departure of the pn spectrum in this 
energy range is present also compared with the RGSl and 
RGS2. Note that this excess is not at the same energy of 




31 32 33 34 35 3{ 

S + H (s"') 

Fig. 2. Hardness ratio (H = 2-10 keV; S = 0.3-2 keV 
counts) versus total intensity (S-l-H, 0.3-10 keV) for the 
two pn observations (squares mark the first observation, 
circles the second). The binning is 512 s. 
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the feature present in the ASCA SIS spectra, which was 
interpreted as due to additional absorption by neutral Ne. 
Smaller differences between the pn and MOS cameras are 
also present at other energies, below 0.4 keV, and up to 
about 1.7 keV. Uncertainties in the calibration of the pn 
Small Window mode below 2 keV are reported in Kirsch et 
al. H2(J04|I . although the use of the latest updated response 
matrices should reduce these differences (Saxton I2004|) . 
Examination of the pn background spectra does not re- 
veal any features at these energies, so we are confident 
that they are not due to improper background subtrac- 
tion. Moreover, both observations show similar shapes for 
the structured residuals. Thus, we restricted the pn en- 
ergy range to 1.7-12 keV, where there are no significant 
differences between the MOSl, M0S2 and pn spectra. 

In summary, we used the following energy ranges: 0.4- 
2 keV for RGSl and RGS2, 0.3-8 keV for the MOSl (be- 
cause of the low statistics at high energy), 0.3-10 keV 
for the M0S2, and 1.7-12 keV for the pn. The RGS 
spectra in both observations do not show evidence for 
edges or emission features. We investigated the 0.3-12 keV 
4U 1850—087 spectra by simultaneously fitting the RGSl, 
RGS2, MOSl, M0S2 and pn spectra of each of the two 
observations individually. Factors were included in the 
spectral fitting to allow for normalization uncertainties 
between the instruments. In the spectral fitting xspec 
version 11.2 was used, and the interstellar abundances of 
Wilms et al. H2()()()|l were used in the photoelectric absorp- 
tion models. 




Energy (keV) 



Fig. 3. Residuals (in units of standard deviations) when 
fitting the 2003 September MOS and pn spectra with a sin- 
gle absorbed power-law model. There is a different struc- 
ture of the residuals in the pn (filled circles), compared 
with MOSl (crosses) and M0S2 (triangles). A similar 
structured excess is present in the 2003 October spectrum. 

LMXB X-ray spectra are generally fit with two compo- 
nent models, a black-body or a multicolor disk-blackbody 
component to account for the low-energy emission (origi- 
nating from the accretion disk or the neutron star surface), 
and a high-energy component which is usually modeled 
with power-law, cut-off power-law, or Comptonized com- 
ponents, to account for the high-energy emission thought 



to be produced in a corona. In order to check if a soft 
component was required by the data we tried first with 
the simplest model consisting of an absorbed power-law 
(photoelectric absorption model phabs in xspec). The 
fits resulted in positive residuals between 0.3-0.7 keV and 
reduced — 1.69 for 1098 degrees of freedom (d.o.f.), 
and = 1.94 for 969 d.o.f. for the first and the second 
observations, respectively. Adding a blackbody improved 
the fits (reduced = 1-51 for 1096 d.o.f, and = 1-78 
for 967 d.o.f.), but structured residuals at low-energies re- 
main. Using a disk-blackbody component instead of the 
blackbody resulted in better fits, with reduced x^=1.49 
for 1096 d.o.f., and x^=1.76 for 967 d.o.f. for the first and 
second observations, respectively. 

The absorption resulting from these fits was always 
higher than that derived from the optical reddening to 
NGC6712. As well as the interstellar absorption in the 
direction of the globular cluster, modeled with PHABS 
with Nh fixed at 1.8x10^^ cm~^, we added another mul- 
tiplicative component to investigate whether an ionized 
absorber could be present (absori model in xSPEC). We 
fixed the Fe abundance of the absorber to the NGC 6712 
value ([Fe/H]=— 0.80 ± 0.2), and linked the photon index 
of the ionizing continuum to that of the power-law com- 
ponent. The fit resulted in an un- ionized cold absorbing 
medium, with a best-fit ionization ^ = L/nB? in the range 
0.5-2 (where L is the ionizing luminosity, n the density of 
the absorbing medium, and R is the distance of the ob- 
scuring material to the source). Thus, we do not consider 
further the presence of any ionized absorber. 

We next tried partial covering (pcfabs model in 
xspec), absorbing both the disk-blackbody and power- 
law continuum components. We included a PHABS com- 
ponent with TVh fixed at 1.8x10^^ cm~^ to account for 
the interstellar absorption. The fit resulted in additional 
absorption in the range 6-8x10^^ cm^^ for the two ob- 
servations, with a covering factor ~95% for reduced 
= 1.37 for 1095 d.o.f., and = 1-40 for 966 d.o.f. for 
the first and the second observations, respectively. The 
disk-blackbody temperatures were 0.6 and 0.3 keV (with 
the innermost radii of the accretion disk rin(cosi)'''^ = 
3.9±0.5 km and 17±6 km) while a = 2.22±0.05 and 
2.35±0.02 (with power-law normalizations of 0.060^Q QQg 
and 0.090±0.003 photons keV"! cm'^ s'^ at 1 keV, for 
the first and the second observations, respectively). 

In order to investigate the Ne/0 abundance ratio, we 
used a variable abundance absorption model (vphabs in 
xspec), with the elemental abundances set to the ISM val- 
ues of Wilms et al. f20DD|) except for those of O, Ne and Fe 
which were fixed to zero. Their absorption effect has been 
replaced with three edges (0-K, Fe-L, Ne-K edges) with 
energies fixed at 0.54, 0.71 and 0.87 keV, and edge depths 
allowed to vary. In this way we could also account for a lo- 
cal iron abundance likely different from the cosmic value. 
We fit the spectra with 3 different two-components mod- 
els for the continuum: (1) a disk-blackbody (diskbb in 
XSPEC, Mitsuda et al. ll984|l and a power-law, (2) a black- 
body and a power-law, and (3) a disk-blackbody and a 
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Table 2. Best-fit continuum parameters when the RGSl, RGS2, pn, MOSl and M0S2 4U 1850-087 spectra were fit 
simultaneously. d66=disk-blackbody, pou'=power-law, ctt=Comptonization model COMPTT in xspec. a is the power- 
law photon index, /cTjn the inner disk temperature, rin(cos z)*^'^ is the inner disk radius for a distance of 8 kpc, i is the 
disk inclination angle and i?bb is the blackbody radius. For the COMPTT model, fcTo is the temperature of the "seed" 
photons, kTc is the electron temperature (fixed at 70 keV) and Tp is the plasma optical depth. 



Model 



(10^ 



(keV) 



rin(cosi)'''^ 
(km) 



(keV) 



^bb 

(km) 



(keV) 



xVdof 



Observation 1 
dbb-|-pow 
bb-|-pow 
dbb-l-ctt 



5.71 



+0.3 
0.2 

6 1+0-2 
"■^-0.2 

4 1+0-1 

^■^-0.1 



2.03 
2.18 



0.06 
+0.04 



0.69; 



-0.02 
-0.02 



0.73: 



2.9 



+0.2 
0.2 



2.5; 



-0.3 
-0.2 



0.48; 



-0.02 
-0.02 



5.1 



[-0.6 
-0.6 



0.23; 



-0.01 
-0.03 



J^-^-0.1 



1.34/1092 
1.37/1092 
1.30/1091 



Observation 2 
dbb-|-pow 
bb-|-pow 
dbb-l-ctt 



5.71 



+0.2 
0.3 

5 Q+°-2 

J.»_0 2 
o fj+0.2 



2.16: 



"+0.06 
-0.06 

2 24+°-"^ 

^•^^-0.04 



0.75 



+0.06 
0.04 



n s/i+010 



1.51 



0.55^ 



2.7 



+ 1.0 
0.7 



n 99+0.03 
'-'■^^-0.02 



1 Q+0.1 



1.43/963 
1.44/963 
1.38/962 



Table 3. Photoelectric absorption results towards 
4U 1850—087 (see TableElfor the continuum emission pa- 
rameters). The edge energies were fixed at 0.54, 0.71 and 
0.87 keV for O, Fe and Ne, respectively. Todgo is the absorp- 
tion depth. A^z is the element column density (in units of 
10^^ cm~^) calculated using the Henke et al. (1993) cross 
sections. TVh is the hydrogen column density implied by 



Nz, in units of 10^ 
of Wilms et al. J^T 



cm , assuming the ISM abundances 



Model 


Edge 


"^cdgc 






Ne/0 


Observation 1 


1Q+0-04 
^■^^-0.05 

2 07+0-12 

^■'J' -0.12 
>J-O*~0.06 




6.oll;^6 

7.4ir4 

9 q+0.3 




dbb-l- 
pow 


Ne K 
O K 

Fe L 


5.2l};i 

36.4+2-1 

o.r7l°;« 


14+0.03 


bb+ 
pow 


Ne K 
O K 

Fe L 


'^•^'-'-0.04 
9 -1 7+0.12 
' -0.12 
n C7+0.06 


5.5li;i 
38.2+2-1 

O.SllH^o 


6.3li;^3 

7 0+0.4 
'•0-0.4 
n+0.3 
■-'•"-0.3 


14+0.03 


dbb-h 
+ctt 


Ne K 
K 

Fe L 


n 1S+0.04 

U-18_o.05 

1 76+0-" 

^- '"-0.13 

46+0-0'' 

"•^"-0.06 


5.0il;i 

31-0+2-5 

66+0-10 


5 7+1-^ 
0- '-1.6 

6 3+0-5 

2.4lg:^ 


16+0.04 
'-'•^^"-o.os 


Observation 2 


n 97+0.05 

'^•^'-o.os 

9 09+0.13 
^•'-'^-0.13 

70+0-07 

"-'-0.07 




8.3lg:^ 

q 7+0.4 
■J- ' -0.4 




dbb-h 
pow 


Ne K 
O K 

Fe L 


7 4+1.4 

' -^-1.4 

40.812:3 

i-ol^:i 


r, ia + 0.04 
U-J^O-0.04 


bb+ 
pow 


Ne K 
O K 

Fe L 


n 90+0.05 
'^•^''-0.05 

2 41+0-12 
7S+0-07 

''-'-0.07 


8.oli;l 

42.4+2-1 

i-ol^:i 


9.2li;^6 

8.7ir4 
n+0.4 

■^■^-0.4 


iq+o.03 

'-'•^^-0.03 


dbb-h 
ctt 


Ne K 
K 

Fe L 


25+"-"-* 
1.90t°;i^ 
0.60t°;?I 


6.9li;i 

QQ K + 2.5 

00 . u 3 9 


7.9+ii 

6 q+o-5 

q 9+0.4 
'-'■^-0.4 


21+0-04 

'-'•^1-0.07 



high-energy Comptonized component (comptt in XSPEC, 
I1994II. This latter model was considered since 



Titarchuk | 

it has been successfully fit to almost all the broad-band 
spectra of the galactic globular cluster LMXBs (Sidoli et 



al. I2()()l|l . Since the BeppoSAX best-fit electron temper- 
ature, kTe, of ~70 keV is well above the XMM-Newton 
upper energy threshold, it was fixed to 70 keV for the fits 
performed here. 

All 3 continuum models fit well and the best-fit con- 
tinuum parameters are given in Table[2| All 3 models give 
similar 0.5-10 keV source luminosities of 1.4 x 10'^^ erg s~^ 
and 1.6 x 10^^ erg for a distance of 8 kpc for the first 
and second observations, respectively. Table 13 gives the 
measured Ne/0 abundance ratios. These are all similar to 
the ISM value of 0.18 of Wilms et al. (2001). The equiv- 
alent column densities due to O, Ne and Fe were then 
estimated by modeling the absorption as occurring from 
3 edges with energies fixed at 0.54, 0.71 and 0.87 keV 
together with a narrow absorption line at an energy of 
0.53 keV (to account for O I ISM absorption). The re- 
sulting edge depths, columns and equivalent hydrogen 
columns are listed in Table 13 and shown in Fig. 01 In Fig. El 
we show the best fit spectra during the two observations. 

As a final test, we tried fitting the spectra with the 
ASCA model, with the Ne and O abundances fixed at 
the ASCA best-fit values (Juett et al. 2001), letting aU 
the other parameters vary in the usual way. We obtained 
unacceptable fits with reduced = 2.0 for 1096 d.o.f., and 

= 2.2 for 967 d.o.f. with structured residuals evident 
below 1 keV (see Fig. El). 



4. Discussion 

We report the results of two XMM-Newton observa- 
tions of 4U 1850—087 performed about 12 days apart. 
The spectra require a soft emission component, which is 
slightly better described by a multi-color disk-blackbody 
than a blackbody. At higher energies, a power-law pro- 
vides a good fit the spectra. The photon index is similar 
to that measured during the ASCA observation (Juett 
et al. I2()()l|l . while it is significantly softer than dur- 
ing the later Chandra observation (Juett & Chakrabarty 
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Fig. 4. The hydrogen column densities (in units of lO'^^ cni~^) impUed by the optical depths of the Ne, O and Fe 
edges, from the first (left) and second observations (right), depending on the different models assumed for the continua 
(Table |3J). Shadowed regions mark the 90% confidence range for N-^ resulting from the contimmm fits (vphabs model, 
see Table El . 



I2()()5|l . This may be due to the fact that in fitting the 
Chandra LETGS spectrum the low-energy absorption was 
fixed to the ISM value towards the globular cluster. We 
note however, that the source luminosities during the 
XMM-Newton and Chandra observations were similar, 
whilst during the ASCA observation the source was al- 
most a factor 2 brighter with a 0.5-10 keV luminosity of 
2.5 X 10^6 erg s'^. 



The total low-energy absorption resulting from the fits 
is similar in both the XMM-Newton observations and is 
4-6.3x10^^ cm~^, depending on the model adopted for 
the continuum. There is evidence for extra-absorption in 
the line of sight, since the best-fit total Nh is always signif- 
icantly higher than the optically derived value in the direc- 
tion of the host globular cluster of (1.8±0.2) x 10^^ cm~^. 
Thus the intrinsic absorption ranges from 2 to 4.5x10^^ 
cm~^, depending on the contimmm model assumed. In 
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Fig. 5. The 0.3-12 keV 4U 1850—087 count spectra (together with the residuals in units of standard deviation) from the 
two observations (left panel shows the first observation, the right panel the second). The continua have been modeled 
with a disk-blackbody and a power-law (see Tables |2] and |31 for the continuum and line parameters, respectively). 
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Fig. 6. The 0.3-12 keV 4U 1850— 087 count spectra (together with the residuals in miits of standard deviations) from 
the two observations (left panel shows the first observation, the right panel the second), fit with the ASCA model with 
the Ne and O abundances fixed at the values reported in Juett et al. ((200l|l . 



the Comptonization model a lower column density is re- 
quired because of the turnover at low energies present in 
the model. The presence of neutral extra-absorption local 
to the source is also confirmed by a good fit when using the 
partial covering fraction absorption model, which indicate 
that the central source is absorbed by a neutral medium 
with a covering factor of ~95% and an intrinsic hydrogen 
column density in the range 6-8x10^^ cm~^ for the two 
observations. 

We adopted a variable absorption model together with 
three edges (0-K, Ne-K and Fe-L) in order to measure the 
column density of the Ne, O and Fe in the line of sight, 
since the ASCA spectrum suggests an excess absorption 
of neutral Ne-rich material local to the source (Juett et 
al. l2()()l|l . Other ultra-compact X-ray binaries (Porb<l hr) 
display over-abundances of neutral Ne from the absorption 
effects in ASCA spectra. Moreover, during XMM-Newton 
and Chandra observations an anomalously high Ne/0 
abundance ratio has been observed in a a number of other 
ultra-compact binaries, indicative of neutral Ne overabun- 



dance (e.g., 4U 0614+091, Pacrcls et al. lMITl 4U 1543-624 
and 2S 0918-549, Juett & Chakrabartv lMI5|l . Thus, it has 
been proposed that the donor stars in some ultra-compact 
binary systems are Ne-rich white dwarfs (e.g., Yungelson 
et al.EnnS Bildsten'2D^). 

The high resolution RGS 4U 1850-087 spectra pre- 
sented here do not show any prominent emission fea- 
tures or other absorption edges, besides those of 0-K, 
Ne-K and Fe-L. Our study of the absorbing Ne and O 
toward 4U 1850— 087 reveals an Ne/0 abundance ratio 
(see Table which is consistent with the ISM value 
of 0.18 (Wilms et al. 12000(1 . This result is in agreement 
with the Chandra observation performed in 2002 (Juett 
& Chakrabarty, 12005(1 but contrary to the earlier ASCA 
measurement (Juett et al. EDOlf) . 

If the measured elemental column densities Nz (see 
Table 13) are converted to equivalent H column densities 
using the Wilms et al. 1(2000(1 ISM abundances, we can 
compare the A^h resulting from the overall shape of the 
XMM-Newton spectra (dashed regions in Fig. 0] include 
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the 90% uncertainties on the A^h resulting from the fit). 
For each continuum model, for both observations, there 
seems to be a discrepancy between the total Nh and the 
equivalent A'^h calculated from the elemental column den- 
sities Nz ■ This could indicate an over-abundance of Ne and 
O and a sub-solar abundance of Fe. This under-abundance 
may be explained by the low Fe abundance of the host 
globular cluster. Alternatively, it is possible that the un- 
certainties on the derived hydrogen column densities could 
be underestimated. Indeed, Paerels et al. (2000) point out 
that the photoelectric cross-sections could have a 30% un- 
certainty which is large enough to account for this discrep- 
ancy. 

There is no evidence for any temporal variability of 
the Ne, O and Fe column densities between the two 
XMM-Newton observations (within 90% uncertainty), al- 
though we note that the column densities are systemati- 
cally higher in the second observation. Also the Ne/0 ra- 
tio, although always compatible within uncertainties with 
the standard ISM value, is on average higher during the 
second observation. The 4U 1850—087 spectrum during 
the two XMM-Newton observations differs in the total 
0.5-10 keV luminosity (see Table |2l with the second ob- 
servation being ~10% more luminous than the first. This 
suggests a possible correlation between the Ne/0 abun- 
dance ratio and the X-ray source luminosity, and a pos- 
sible explanation for the different Ne/0 ratios observed 
with XMM-Newton, Chandra, and ASCA. Note that dur- 
ing the ASCA observation, where evidence for a strong Ne 
overabundance was reported, the source luminosity was al- 
most twice that during the XMM-Newton and Chandra 
observations. 

Juett & Chakrabarty H2003|l proposed that ionization 
could play a role in the variable abundance ratios. We 
suggest another possible mechanism which could help in 
the understanding of the local extra-absorption and its 
metal abundance. Maccarone et al. (2004) studied the 
irradiation-induced stellar winds in X-ray binaries in or- 
der to explain different X-ray spectra from LMXBs lo- 
cated in globular clusters with different metallicities. An 
evaporative wind can be produced even in LMXBs with 
degenerate companions, where part of the radiation pro- 
duced from the central source illuminates the donor star 
(e.g., Ruderman et al. I1989|l . This wind could contribute 
to the observed column density toward the X-ray source, 
leading to an intrinsic column density of ~6xl0^^ cm~^ 
(see eq. (7) in Maccarone et al. OITOji . During the XMM- 
Newton observations we found a comparable amount of 
extra-absorption towards 4U 1850— 087. We suggest that 
a wind evaporated from the degenerate companion could 
be responsible for the intrinsic absorption observed. Using 
the 4U 1850—087 observed parameters, and eq. (7) of 
Maccarone et al. ()2()04|l . we derive a wind velocity of 
~5xl0'' cm s^^, which may be confined to the binary 
(the escape velocity is ~10® cm s""'^). A higher source 
luminosity would translate into a larger contribution by 
the wind from the degenerate donor, which is likely to be 
rich in Ne and O. We suggest that this mechanism could 



contribute to the different abundance ratio observed from 
4U 1850-087 with XMM-Newton, Chandra, and ASCA. 
This could possibly help in explaining why during higher 
source luminosity intervals, higher Ne/0 abundance ratios 
are observed. We note that another ultra-compact X-ray 
binary, 4U 1543-624, displays a variable Ne/0 abundance 
ratio. Juett & Chakrabarty (2003) measured an Ne/0 
abundance ratio of 1.5±0.3 with Chandra, and 0.54±0.03 
with XMM-Newton. The higher Ne/0 abundance ratio 
was observed when the source was more luminous, as ap- 
pears to be the case with 4U 1850—087. 

Theoretical models for the formation of white dwarfs 
predict a Ne/0 abundance ratio in the range 0.2-0.4 (e.g., 
Deloye & Bildsten IM)^ Segretain et al. IBMl Gutierrez 
et al. rr99f?jl- which is low compared with the Ne/0 ratios 
observed with ASCA in 4U 1850-087 (or with Chandra in 
4U 1543-624). On the other hand, Yungelson et al. 12001 
studying the formation of Ne-enriched donors in ultracom- 
pact X-ray binaries, point out that the abundance of neon 
in the nucleus of the dwarf may be well underestimated by 
a factor of 3 (Isern et al. 11991(1 . which makes the theoret- 
ically predicted Ne/0 ratios agree better with the higher 
observed values. 
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